Real-time dynamics of the O(N) model in 1+1 dimensions 
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We study the non-equilibrium dynamics of the O(N) model in classical and quantum field theory in 
1+1 dimensions, for N > 1. We compare numerical results obtained using the Hartree approximation 
and two next to leading order approximations, the bare vertex approximation and the 2PI-1/N 
expansion. The later approximations differ through terms of order g 2 , where g is the scaled coupling 
constant, g = X/N. In this paper we investigate the statement regarding the convergence with 
respect to g. We find that the differences between these two approximation schemes diminish for 
larger values of N, when A is fixed. 
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I. INTRODUCTION 

Initial value problems in quantum field theory repre- 
sent a key step in the pursuit of first principle under- 
standing of physics of the early universe and dynamics 
of phase transitions and particle production following a 
relativistic heavy-ion collision. In order to achieve quan- 
titative comparisons with experimental data, it is abso- 
lutely necessary to go beyond existing mean field theo- 
retical models, and implement approximations which ac- 
count for higher order corrections. Thus we will not only 
have a reliable picture of the phase diagram of quan- 
tum field theories, but the model will include important 
physics such as a mechanism of allowing an out of equilib- 
rium system to be driven back to equilibrium, mechanism 
which is absent in leading order. Much progress has been 
achieved in recent years (see Ref.r:BS03, and references 
therein.) 

We have recently 0] presented a first study of the dy- 
namics of a single component "explicitly broken sym- 
metry" A <p A field theory in 1+1 dimensions, using a 
Schwinger-Dyson (SD) equation truncation scheme based 
on ignoring vertex corrections we call the bare vertex ap- 
proximation (BVA) 0| ■ We have compared the BVA pre- 
diction for the time evolution of the system with results 
obtained using two other approximations: the Hartree 
approximation which represents a leading order (LO) re- 
sult in our three-point vertex function truncation scheme, 
and the 2PI-1/N expansion Q, which is a related next 
to leading order (NLO) approximation scheme. The 2PI- 
1 /N expansion can be formally obtained from the BVA by 
dropping a graph from the two-particle irreducible (2PI) 
effective action 00,0|. 

Two important aspects in connection with this previ- 
ous study deserve attention: Firstly, we have observed 
that the BVA approximation of the 0(1) model in 1+1 
dimensions offers a unique toy problem which allows for 
studying the dynamics of phase transitions. It is well 
known that in 1+1 dimensions there is no phase transi- 
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tion in the 0(1) model except at zero temperature 0. 
Nevertheless, we used this model to demonstrate some of 
the features expected to be true in 3+1 dimensions such 
as the restoration of symmetry breakdown at high tem- 
peratures and equilibration of correlation functions. For 
this model, we have found evidence that the order of the 
phase transition predicted by the static phase diagram in 
the Hartree approximation, is relaxed from a first order 
to second order phase transition in the BVA. Moreover, 
unlike in the case of the Hartree approximation where 
the order parameter oscillates about the Hartree mini- 
mum without equilibration, the BVA results track the 
Hartree curve, except with damping: the average fields 
in the BVA equilibrate for arbitrary initial conditions. 
The 2PI-1/N expansion goes to a zero value of the or- 
der parameter, in agreement with the exact result of no 
phase transition for finite temperature. 

Secondly, we have shown that the disagreement be- 
tween the 2PI-1/N and BVA predictions for the order 
parameter is more pronounced at larger values of the 
coupling constant A, while predictions for the effective 
mass are similar for both approximations, in agreement 
with expectations based on our experience with the clas- 
sical limit of these approximations 0]- We remind the 
reader that in the classical field theory case in 1+1 di- 
mensions we were able to compare these approximations 
with exact Monte Carlo simulations 0| . Note also that in 
the classical case the quality of the agreement between 
approximation and exact results is not sensitive to the 
coupling constant value, since the coupling can be scaled 
out in the classical limit. This scaling is not possible 
in the quantum case. Nevertheless, direct comparisons 
with exact results are currently possible only in the clas- 
sical case, and the differences between the BVA and the 
exact results were said to be due to the fact that three- 
point vertex corrections become important at interme- 
diate time scales, and they affect the evolution of the 
two-point function more than the time evolution of the 
order parameter. 

The study of the 0(1) model in 1+1 dimensions car- 
ried out in was both exciting and worrisome: excit- 
ing because the model provided a testing ground where 
features of the theory expected to be true in 3+1 di- 
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mcnsions could be studied; worrisome because the huge 
discrepancy between the 2PI-1/N and the BVA results 
made it look hopeless to think that one can accurately 
predict the time evolution of the system. It is imper- 
ative to achieve convergence with the truncation order 
of the SD hierarchy of equations, in order to allow for 
reliable predictions of the dynamics of the system, and 
the results for the 0(1) model in 1+1 dimensions imply 
that the difference between the 2PI-1/N and BVA ap- 
proximations is considerable. Since the formal difference 
is related to the additional graph in the 2PI effective ac- 
tion considered in the BVA, the importance of this graph 
appeared to be considerable. Since the physical system 
we are interested in is in fact the 0(4) linear sigma model 
in 3+1 dimensions, one can only hope that by increasing 
the number of dimensions and/or the number of fields, 
the agreement between these methods will get better. In 
this paper wepresent an immediate extension of the work 
discussed in U , by generalizing the calculations from the 
0(1) to the 0(N) model, in 1+1 dimensions. 



II. BVA EQUATIONS 

In this section we briefly review the basics of the BVA 
formalism (for full details, see Ref. @): We start with 
the Lagrangian for the 0(N) symmetry model 
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(i) 



Here the Einstein summation convention for repeated in- 
dices is implied. We introduce a composite field x( x ) by 
adding to the Lagrangian a term: 



This gives a Lagrangian of the form 



(2) 
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A 2A ' 

which leads to the classical equations of motion 

□ +X(a:)] =0, 
and the constraint ( "gap" ) equation 



A 



(3) 



(4) 



(5) 



It is convenient now to introduce the super-fields and 
super-currents, (f> a (x) = {x(x),<f) 1 (x),^ 2 (x),. . .,<f> N (x)}, 
and j a (x) = {J(x),j 1 (x),j 2 (x),...,j N (x)}, with a = 



0,1, ... ,N. Then, the BVA equations can be obtained 
from the effective action 

T[0,G] =5 cl [0] + ^Trln[G- 1 ] + ^Tr[G o - 1 G]+r 2 [G], (6) 

where T2[G] is the Cornwall- Jackiw-Tomboulis gener- 
ating functional of the two-particle irreducible (2-PI) 
graphs. Here G denotes the Green function matrix 

' X) - S J0 (x')-{K t (x,x') G^x,x')) ' (7) 
and iS c i[0] is the classical action in Minkowski space 

S44>] = J d 2 x{~^Mx)[a + x (x)}Mx) 



X 2 {x) . M 



2A 



+ T*<*>} 



(») 



In the BVA approximation, we keep in T 2 [G] only the 
terms 



r 2 [G] 



dx dy 



Gij(x,y)Gji(y,x)D(x,y) (9) 
+ 2K i (x,y)K j (x,y)G ij (x,y) . 



This is equivalent to a truncation of the SD tower of 
equations, where one approximates the exact three-point 
vertex function equation 



r a /3 7 [(/>] (Xl ,X 2 ,X 3 ) =f a 0-y Sc {xi ,x 2 )Sc (xi , x 3 ) 
SE a/3 [(j)\(xi,X 2 ) 



H-t{x 3 ) 



(10) 



Here f ij0 = /, 



f. 



ij , and / is zero otherwise, 
and T, a p[(j)}(x,x') denote the elements of the self-energy 
matrix. Then, the Hartrce approximation corresponds 
to the case when one ignores the three-point vertex cor- 
rections altogether, while in the BVA one keeps only the 
contact term. Finally, the 2PI-1/N approximation is ob- 
tained from the BVA approximation by dropping the sec- 
ond term in r 2 [G], see Eq.JTU}. 
The following equations of motion 

n + x (x)]<k(x)+K i (x,x)/i = 0, (11) 



X(x) 



^2[4>i{ x ) +Gu{x,x)/i 



(12) 



must be solved sclf-consistently with the Green function 
equations 

G~l{x, x') = G-l p {x, x 1 ) + Z af3 (x, x 1 ) , (13) 
which are subject to the condition 

dx" G-^x, x") G M {x",x') = 6 ay 5 c (x, x') . (14) 
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Here g = X/N, and this indicates that the expansion pa- 
rameter in this formalism is g, rather then the usual 1/N 
corresponding to the familiar large-N expansion. The in- 
tegrals and delta functions 8c{x,x') are defined on the 
closed time path (CTP) contour, which in corp orates the 
initial value boundary condition 0, 0, Mj. The 



elements of the Green function G Q ^{x, x') arc 



G 



00 



r -l 



— 6c(x,x') , 
9 



(15) 



G (x, x') 

K oj( x ' x ') = ^ji^^ci^x') , 
G^ = K^(x,x') = K^(x\x). 

For the BVA, the self-energies £^J A (x, x') are given as 

Il(x,x') = ^G ij {x,x')G ij (x,x'), (16) 

£li(x, x') = i Kj(x, x') Gji(x, x') , 
Cli(x, x) = i Gijix, x) Kj(x, x) , 

Tiij(x,x') =i Gij(x,x ) D(x,x ) + Ki(x,x ) Kj(x,x ) 

The equations for the 2PI-1/N approximation can be for- 
mally obtained by setting Qi(x,x') = Qj(x',x) = 0, and 
dropping the KiKj term in the definition of x'). 

The BVA is an energy-conserving approximation, 
where the average energy is given as 



E 



' A 



<X> - ^(X 2 ) + \{{[dtcp} 2 ) + {[d x 4>] 2 ) + <X</> 2 )} , 

(17) 

with the expectation values 

(X 2 ) = X 2 (x)+D(x,x)/i, 
(0?) =$(x) + Gu(x,x)/i, 
(X<t> 2 ) =x{x)Y J ^ 2 ) + 2Mx)Ki(x,x)/i 

i 

+ \ [ dx'\u(x,x')D(x',x) + Cl i (x > x')K i (x',x) 
3 Jc L 

+ Sli(x,x')Ki(x',x) + T>ij(x,x')Gji(x',x) . 



(18) 
(19) 
(20) 



Since the energy calculation involves only second order 
derivatives of the effective action, it results that any trun- 
cation scheme done at the level of the three-point vertex 
function or beyond, and which produces a self-energy 
matrix obeying the appropriate symmetries, will be an 
energy-conserving approximation. 

In order to investigate the convergence of the proposed 
approximations with respect to g, we consider two initial 
scenarios which, from the point of view of computational 
storage requirements, differ only minimally from the 0(1) 
calculations we have previously performed. We summa- 
rize next the explicit form of the equations we have to 
solve, and propose a practical way of obtaining the nu- 
merical solution. 



A. One-field scenario 

Firstly, in what we call the "one-field scenario" , we 
study the case when only one of the initial fields, denoted 
4>\, is nonzero, i.e. 2 (O) = • • • = <^jv(0) = 0, and <fo(0) = 
• • • = 0jv(O) = 0. In this case, we can show that the 
Green function matrix remains diagonal at all times, and 
all diagonal matrix elements other than Gn are identical. 
Therefore, in the one-field scenario the Green function 
and self-energy matrices can be written as 



(21) 
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(22) 



(23) 



We find it convenient to introduce the additional Green 
functions 

D^ l {x,x') = Dq 1 {x,x') +U(x,x') , (24a) 
Gi 1 (x,x') = Gq X {x,x') + Sii(.x,x') , (24b) 
Gz 1 (x,x')=G - 1 {x 1 x') + J: 22 (x,x'), (24c) 
~o(x,x') = Ka x {x,x') + Q.{x,x') . (24d) 

With these notations, the Green function equations are 

l = D^D + EoK, (25) 
= D^ 1 K + E G 11 , 

= EoD + G^K , 

1 = S K + G\ 1 Gn , 

'22 • 



1 — G 2 1 G2 



The numerical solution is achieved by iterating the fol- 
lowing system of equations 



K = -D 2 E Gn , 
K = - Gi S D , 

D = D 2 + D 2 [s Gi S 
Gn = Gi + Gi D 2 So 
G22 = G2 . 



(26) 



D, 
Gn 
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B. N-copies scenario 



III. RESULTS 



Secondly, we study the case when initially we have N 
identical copies of the field, i.e. <j>i(0) = (j>o an( A 4>i(0) = 
7To, for i = 1 . . . N. In this case we can still choose the 
Green function to be diagonal at t = 0, but this property 
is lost as we propagate the equations of motion. How- 
ever, the Green function matrix in the "N-copies sce- 
nario" remains particularly simple, as the diagonal and 
off-diagonal matrix elements are separately identical. 

Once again, we start with the Green function and self- 
energy matrices, which in this case have the form 
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(27) 



(28) 



(29) 



Correspondingly, the Green functions equations we have 
to solve are 



I = D? 1 D + NZ K , 

Gn + (N-l)Gi2 



ZA7 1 K + Eo 



J 2 

0-- 2 

= Z D 

1 = S A 



(30) 



K, 



Gr i + (N-l)Ei2 

Gr 1 Gn + (N-l)Si 2 Gi2, 
= S K + G- 1 G12 + E 12 Gn + (N - 2)Si2 G 12 



In our strategy, it is convenient to introduce the addi- 
tional set of Green functions 



G- 1 =G^ -E D 2 E , (31a) 
Gr 1 - E12 , (31b) 
Gr 1 + (N- l)Eia . (31c) 



G7 1 



Then, we obtain the desired Green functions by solving 
the equations 
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We study the O(N) model, in 1+1 dimensions. We 
set A = 7.3 since this particular value of A was used 
in previous studies of the dynamics of disoriented chiral 
condensates (DCC) in 3 + 1 dimensions in the leading 
order in large-N approximation |l4l Il5j . We remind the 
reader that g = X/N, and this is the actual expansion 
parameter in this approximation. 

We assume here that the initial state is described by 
a Gaussian density matrix peaked around some non-zero 
value of {(f>(0)), and characterized by a single particle 
Bose-Einstein distribution function at a given tempera- 
ture. The details of this choice of initial conditions has 
been discussed in detail in Ref. Q]. In this study, we 
choose a very low initial temperature of To = 0.1 in or- 
der to emphasize quantum effects in the dynamics. We 
set the initial values of the non-vanishing initial fields to 
0o = 4>(0)/N = 1, and tt = 0(0) = 0. 

The numerical procedure for solving the BVA equa- 
tions is described in detail in Refs. (16tfl7j| . The unknown 
functions are expanded out in terms of Chcbyshcv poly- 
nomials, on a nonuniform grid. The algorithm follows a 
multi-step approach, and the algorithm possess spectral 
convergence. We achieve convergence of the numerical re- 
sults with a relatively small number of grid points, only 
32 and 128 points, for the time step and the momentum 
domain discretization, respectively. The self-energy ma- 
trix elements are calculated using standard fast-Fourier 
transform algorithms on a uniform 1024 points grid. We 
employ a cubic-spline interpolation technique to perform 
the necessary transformation between the two momen- 
tum grids. The mass renormalization necessary for the 
quantum field theory calculations has been discussed pre- 
viously 0]. A momentum cutoff, A = 37r, was chosen for 
the purpose of the present calculation. 

We study the time evolution of the one-point func- 
tions: the order parameter {(f>(t)). and the auxiliary field 
(x(t))- We compare the Hartree, 2PI-1/N expansion and 
the BVA approximation. Results are presented for both 
the classical (CFT) and the quantum field theory (QFT) 
case. In keeping with our earlier work, the potential is 
a double-well potential for the quantum case, but single- 
well in the classical case. 

We start by reviewing the results for the 0(1) model. 
In Fig. n we present the time evolution of {<fi(t)) and 
(x(t)) for CFT and QFT. 

Figures [2] and [3] present the time evolution of the one- 
point functions for the one-filed scenario corresponding 
to the 0(2) and 0(4) models, for the CFT and QFT case, 
respectively. Similar results are shown in Figs. ^ and [S] 
for the N-copies scenario. 

One striking feature in reviewing these results, is the 
fact that the NLO corrections are more important in the 
quantum than in the classical case. The Hartree result 
represents a poor approximation in QFT, while in CFT 
it is arguably doing pretty well. This trend is still appar- 
ent in the N = 10 results (not shown). For N = 10, the 
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FIG. 1: N=l; legend: BVA (solid) , 2PI-1N (dashed), and 
Hartrcc (dotted) 



FIG. 4: N-copies scenario: classical field theory. 
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FIG. 2: One-field scenario: classical field theory. 
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FIG. 3: One-field scenario: quantum field theory. 



differences between 2PI-1/N and BVA for the QFT case 




FIG. 5: N-copies scenario: quantum field theory. 



have essentially disappeared, but Hartree is still consid- 
erably different from NLO in the QFT case. 



IV. CONCLUSIONS 

We have shown here first calculations of the O(N) 
model in classical and quantum field theory of the BVA 
and 2PI-1/N for N > 1 in 1+1 dimensions. For JV = 1, 
earlier calculations had shown substantial differences be- 
tween these two approximations for the quantum field 
theory case. Our results indicate that these differences 
diminish for larger values of N. The BVA and 2PI-1/N 
approximations differ through contributions which are 
proportional to g 2 . Since g = X/N, and we chose to keep 
the coupling constant A fixed, in order to approach the 
limit of the physically interesting realization of this model 
(see Refs. [lj,[l5|), this trend was expected, and we have 
in fact observed it earlier in our previous studies of the 
quantum mechanical version of this model Q ■ What is 
gratifying is the fact that already for N = 4 the two ap- 
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proximations agree rather well. This gives us hope that 
our future calculations of the 0(4) linear model in 3+1 
dimensions will provide reliable predictions, converged 
with respect to the approximation order. Of course, fur- 
ther numerical evidence will be obtained by carrying out 
a similar study by carrying out simulations in 2+1 and 
3+1 dimensions. An additional study will be pursued 
in order to provide a better treatment of the three-point 
vertex, a study which is clearly necessary in order to re- 
duce the differences between the exact and BVA results 
in 1+1 dimensional classical field theory. It is important 
to remember that no exact calculations are available for 
QFT, unlike the CFT case where lattice simulations have 
been performed. Therefore, the fact that the BVA and 
2PI-1/N agree, does not necessarily imply that these ap- 



proximations also agree with the exact result, and nu- 
merical evidence supporting the fact that these further 
corrections do not affect the results in the 0(N) model 
for N > 1, is critical to assert the credibility of the results 
obtained in this approach. 
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